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ABSTRACT 

Aims. The test-field method for computing turbulent transport coefficients from simulations of hydromagnetic flows is 
extended to the regime with a magnetohydrodynamic (MHD) background. 

Methods. A generalized set of test equations is derived using both the induction equation and a modified momentum 
equation. By employing an additional set of auxiliary equations, we derive linear equations describing the response 
of the system to a set of prescribed test fields. Purely magnetic and MHD backgrounds are emulated by applying an 
electromotive force in the induction equation analogously to the ponderomotive force in the momentum equation. Both 
forces are chosen to have Roberts flow-like geometry. 

Results. Examples with an MHD background are studied where the previously used quasi-kinematic test-field method 
breaks down. In cases with homogeneous mean fields it is shown that the generalized test-field method produces the 
same results as the imposed-field method, where the field-aligned component of the actual electromotive force from 
the simulation is used. Furthermore, results for the turbulent diffusivity tensor are given, which are inaccessible to 
the imposed-field method. For MHD backgrounds, new mean-field effects are found that depend on the occurrence of 
cross-correlations between magnetic and velocity fluctuations. For strong imposed fields, a is found to be quenched 
proportional to the fourth power of the field strength, regardless of the type of background studied. 

Key words, magnetohydrodynamics (MHD) - turbulence - Sun: magnetic fields - stars: magnetic fields 



1. Introduction 

Astrophysical bodies such as stars with outer convective 
envelopes, accretion discs, and galaxies tend to be mag- 
netized. In all those cases the magnetic field varies on a 
broad spectrum of scales. On small scales the magnetic 
field might well be the result of scrambling an existing 
large-scale field by a small-scale flow. However, at large 
magnetic Reynolds numbers, i.e. when advection dominates 
over magnetic diffusion, another source of small-scale fields 
is small-scale dynamo action (Kazantsev 1968). This pro- 
cess is now fairly well understood and confirmed by nu- 
merous simulations (Cho & Vishniac 2000; Schekochihin et 
al. 2002, 2004; Haugen et al. 2003, 2004); for a review see 
Brandenburg & Subramanian (2005). Especially in the con- 
text of magnetic fields of galaxies, the occurrence of small- 
scale dynamos may be important for providing a strong 
field on short time scales (10 7 yr), which may then act as a 
seed for the large-scale dynamo (Beck et al. 1994). 

In contemporary galaxies the magnetic fields on small 
and large length scales are comparable (Beck et al. 1996), 
but in stars this is less clear. On the solar surface the so- 
lar magnetic field shows significant amounts of small-scale 
fields (Solanki et al. 2006). The possibility of generating 
such magnetic fields locally in the upper layers of the con- 
vection zone by a small-scale dynamo is sometimes referred 
to as surface dynamo (Cattaneo 1999; Emonet & Cattaneo 
2001; Vogler & Schiissler 2007). On the other hand, simu- 
lations of stratified convection with shear show that small- 
scale dynamo action is a prevalent feature of the kinematic 



regime, but becomes less important when the field is strong 
and has saturated (Brandenburg 2005a; Kapyla et al. 2008). 

An important question is then how the primary pres- 
ence of small-scale magnetic fields affects the generation 
of large-scale fields if these are the result of a dynamo pro- 
cess that produces magnetic fields on scales large compared 
with those of the energy-carrying eddies of the underlying 
and in general turbulent flow (Parker 1979) via an insta- 
bility. A commonly used tool for studying these large-scale 
dynamos is mean-field electrodynamics, where correlations 
of small-scale magnetic and velocity fields are expressed in 
terms of the mean magnetic field and the mean velocity us- 
ing corresponding turbulent transport coefficients or their 
associated integral kernels (Moffatt 1978; Krause & Radler 
1980). The determination of these coefficients (e.g., a ef- 
fect and turbulent diffusivity) is the central task of mean- 
field dynamo theory. This can be performed analytically, 
but usually only via approximations which are hardly justi- 
fied in realistic astrophysical situations where the magnetic 
Reynolds numbers, ReM, are large. 

Obtaining turbulent transport coefficients from direct 
numerical simulations (DNS) offers a more sustainable al- 
ternative as it avoids the restricting approximations and 
uncertainties of analytic approaches. Moreover, no assump- 
tions concerning correlation properties of the turbulence 
need to be made, because a direct "measurement" of those 
properties is performed in a physically consistent situation 
emulated by the DNS. The simplest way to accomplish such 
a measurement is to include, in the DNS, an imposed large- 
scale (typically uniform) magnetic field whose influence on 
the fluctuations of magnetic field and velocity is utilized to 
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infer some of the full set of transport coefficients. We refer 
to this technique as the imposed-field method. 

A more universal tool is offered by the test-field method 
(Schrinner et al. 2005, 2007), which allows the determi- 
nation of all wanted transport coefficients from a single 
DNS. For this purpose the fluctuating velocity is taken 
from the DNS and inserted into a properly tailored set of 
test equations. Their solutions, the test solutions represent 
fluctuating magnetic fields as responses to the interaction 
of the fluctuating velocity with a set of properly chosen 
mean fields. These mean fields will be called test fields. For 
distinction from the test equations, which are in general 
also solved by direct numerical simulation, we will refer to 
the original DNS as the main run. This method has been 
successfully applied to homogeneous turbulence with helic- 
ity (Sur et al. 2008, Brandenburg et al. 2008a), with shear 
and no helicity (Brandenburg et al. 2008b), and with both 
(Mitra et al. 2009). 

A crucial requirement on any test-field method is the 
independence of the resulting transport coefficients on the 
strength and geometry of the test fields. This is immediately 
clear in the kinematic situation, i.e., if there is no back- 
reaction of the mean magnetic field on the flow. Indeed, for 
given magnetic boundary conditions and a given value for 
the magnetic diffusivity, the transport coefficients must not 
reflect anything else than correlation properties of the ve- 
locity field which are completely determined by the hydro- 
dynamics alone. For this to be guaranteed the test equa- 
tions have to be linear and the test solutions have to be 
linear and homogeneous in the test fields. 

Beyond the kinematic situation the same requirement 
still holds, although the flow is now modified by a mean 
magnetic field occurring in the main run. (Whether it is 
maintained by external sources or generated by a dynamo 
process does not matter in this context.) Consequently, 
the transport coefficients are now functions of this mean 
field. It is no longer so obvious that under these circum- 
stances a test-field method with the aforementioned lin- 
earity and homogeneity properties can be established at 
all. Nevertheless, it turned out that the test-field method 
developed for the kinematic situation gives consistent re- 
sults even in the nonlinear case without any modification 
(Brandenburg et al. 2008c). This method, which we will 
refer to as "quasi-kinematic" is, however, restricted to situ- 
ations in which the magnetic fluctuations are solely a con- 
sequence of the mean magnetic field. (That is, the primary 
or background turbulence is purely hydrodynamic.) 

The power of the quasi-kinematic method was demon- 
strated based on a simulation of an a 2 dynamo where the 
main run has reached saturation with mean magnetic fields 
being Beltrami fields (Brandenburg et al. 2008c). Magnetic 
and fluid Reynolds numbers up to 600 were taken into 
account, so in some of the high ReM runs there was cer- 
tainly small-scale dynamo action, that is, a primary mag- 
netic turbulence bo should be expected. Nevertheless, the 
quasi-kinematic method was found to work reliably even 
for strongly saturated dynamo fields. This was revealed by 
verifying that the analytically solvable mean-field dynamo 
model employing the values of a and turbulent diffusivity 
as derived from the saturated state of the main run indeed 
yielded a vanishing growth rate. Very likely the small-scale 
dynamo had saturated on a low level so the contribution 
to the mean electromotive force, which was not taken into 



account by the quasi-kinematic method, could not create a 
marked error. 

Limitations of the quasi-kinematic test-field method 
were recently pointed out by Courvoisier et al. (2010) and 
will be commented upon in more detail in the discussion 
section. Indeed, the purpose of our work is to propose 
a generalized test-field method that allows for the pres- 
ence of magnetic fluctuations in the background turbulence. 
Moreover, its validity range should cover dynamically effec- 
tive mean fields, that is, situations in which the velocity and 
magnetic field fluctuations are significantly affected by the 
mean field. 

With a view to this generalization we will first recall the 
mathematical justification of the quasi-kinematic method 
and indicate the reason for its limited applicability (Sect. 
2). In Sect. 3 the foundation of the generalized method will 
be laid down in the context of a relevant set of model equa- 
tions. In Sect. 4 results will be presented for various combi- 
nations of hydrodynamic and magnetic backgrounds having 
Roberts-flow geometry. The astrophysical relevance of our 
results and the connection with the work of Courvoisier et 
al. (2010) will be discussed in Sect. 5. 

2. Justification of the quasi-kinematic test-field 
method and its limitation 

In the following we split any relevant physical quantity F 
into mean and fluctuating parts, F and /. No specific av- 
eraging procedure will be adopted at this point; we merely 
assume the Reynolds rules to be obeyed. Furthermore, we 
split the fluctuations of magnetic field and velocity, b and 
it, into parts existing already in the absence of a mean 
magnetic field, bo and Uq (together they form the back- 
ground turbulence), and parts vanishing with B, denoted 
by bg and Ug. We may split the mean electromotive force 
£ = u x b likewise and get 

£ = £ + £ B (1) 

with 



£ = Mo x b , (2) 
£g = u x bg + Ug x b + Ug x bg . (3) 

Note that we do not restrict bg and Ug, and hence also not 
£g to a certain order in B. 

In the present Section we assume that the background 
turbulence is purely hydrodynamic, that is, bo = and 
hence b = bg or, in other words, the magnetic fluctuations b 
are entirely a consequence of the interaction of the velocity 
fluctuations u with the mean field B. 

In a homogeneous medium, the induction equations for 
the total and mean magnetic fields as well as for the mag- 



netic fluctuations read 
FiTi 

-^j- — 7]V 2 B + curl (U x B), (4) 

aJ5 

— = ?yV 2 B + curl (C7 x B + £), (5) 
at 

— = ?/V 2 b + curl (U x b + ux B + £'), (6) 



with £' = uxb — u x b. The solution of the linear equation 
© for the fluctuations b, considered as a functional of u. 
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U and B, is linear and homogeneous in the latter and the 
same is true for 



£g = S~uxbg~uxb. (7) 

If the velocity is influenced by the mean field, that is, if 
u and U depend on B, £ considered as_a functional of B, 
£{B}, is of course nonlinear. However, £, again considered 
as a functional of it, U and B, £{u,U , B}, is still linear 
in B. 

The major task of mean-field theory consists now just in 
establishing a linear and homogeneous functional relating 
£ and B. Making the ansatz 

£ = aB rjVB, (8) 

with VB being the gradient tensor of the mean magnetic 
field, this task coincides with determining the tensors a and 
77, which are of course functionals of u and U. Because of 
linearity and homogeneity we are entitled to employ for this 
purpose various arbitrary vector fields B T (i.e., test fields) 
in place of B in ©, keeping the velocity of course fixed. 
Each specific assignment of B T yields a corresponding 6 T 
and via that an ~£ T and it establishes (up to) three linear 
equations for the wanted components of a. and rj. Hence, 
choosing the number of test fields in accordance with the 
number of the tensor components, and specifying the ge- 
ometry of the test fields "sufficiently independent" from 
each other, the components of a and 77 can be determined 
uniquely. In doing so, the amplitude of the test fields clearly 
drops out (Brandenburg et al. 2008b). 

Is the result affected by the geometry of the test fields? 
An ansatz like (J8J) is in general not exhaustive, but re- 
stricted in its validity to a certain class of mean fields, here 
strictly speaking to stationary fields which change at most 
linearly in space. Consequently, the geometry of the test 
fields is without relevance just as long as they are taken 
from the class for which the £ ansatz is valid, but not for 
other choices. 

For many applications it will be useful to generalize the 
test-field method such that all employed test-fields are har- 
monic functions of position, defined by one and the same 
wavevector k. The turbulent transport coefficients can then 
be obtained as functions of fc and have to be identified 
with the Fourier transforms of integral kernels which de- 
fine the in general non-local relationship between £ and 
B (Brandenburg et al. 2008a). Quite analogous, the in 
general also non-instantaneous relationship between these 
quantities can be recovered by using harmonic functions 
of time for the test-fields. The coeffcients, then depend- 
ing on the angular frequency u>, represent again Fourier 
transforms of the corresponding integral kernels (Hubbard 
& Brandenburg 2009). 

If u and U are taken from a series of main runs with 
a dynamically effective mean field of, say, fixed geometry, 
but from run to run differing strength B, a and r] can be 
obtained as functions of B. Thus, it is possible to determine 
the quenched dynamo coefficients basically in the same way 
as in the kinematic case, albeit at the cost of multiple nu- 
meric work. 

Let us now relax the above assumption on the back- 
ground turbulence and admit additionally a primary mag- 
netic turbulence bo- For the sake of simplicity we will not 
deal here with £$, so let us assume that it vanishes. In the 



representation ((3J) of £ we now combine the first and last 
terms using u = Uq + Ug and obtain 



£ = u x bg + u s x 6 , (9) 



differing from ((7]) by the additional contribution, Ug x bo- 
The quasi-kinematic method necessarily fails here even 
when modifying © appropriately to form an equation for 
bg as only the term u x bg is provided. Obviously, a valid 
scheme must treat also Ug in a test-field manner similar 
to bg. The equation to be employed for Ug has of course 
to rely upon the momentum equation. Due to its intrinsic 
nonlinearity, however, a major challenge consists then in en- 
suring the linearity and homogeneity of the test solutions 
in the test fields. 

3. A model problem 

3.1. Motivation 

We commence our study with a model problem that is sim- 
pler than the complete MHD setup, but nevertheless shares 
with it the same mathematical complications. We drop 
the advection and pressure terms and adopt for the dif- 
fusion operator simply the Laplacian (and a homogeneous 
medium). Thus, there is no constraint on the velocity from 
a continuity equation and an equation of state. However, 
as in the full problem, we allow the magnetic field to exert 
a Lorentz force on the fluid velocity. We also allow for the 
presence of an imposed uniform magnetic field B lmp to en- 
able a determination of the a effect independently from the 
test-field method. The magnetic field is hence represented 
as B — Simp + V x A where A is the vector potential 
of its non-uniform part. The resulting modified momentum 
equation for the velocity U and the (original) induction 
equation read then 

8U 

— = JxB + F K + ^V 2 C7, (10) 
8A 

— = U x B + F M + V V 2 A, (11) 

where we have included the possibility of both kinetic and 
magnetic forcing terms, Fk and -Fm, respectively. (In this 
paper we use hydrodynamic and kinetic forcing synony- 
mously.) We have adopted a system of units in which B has 
the dimension of velocity namely B := B / y/JTp, p =const. 
Defining the current density as J = V x B, it has the 
unit of inverse time. The electric field has then the unit of 
squared velocity. Furthermore, v is the kinematic viscosity 
and rj the magnetic diffusivity. 

As will become clear, the major difficulty in defining a 
test-field method for MHD or purely magnetic background 
turbulence is caused by bilinear (or quadratic) terms like 
J x B and U x B. Hence, taking the U ■ VC/ nonlinearity 
into account would not offer any new aspect, but would blur 
the essence of the derivation and the clear analogy in the 
treatment of the former two nonlinearities. The treatment 
of the advective term follows the same pattern, as is demon- 
strated in Appendix [XJ Also, it should be pointed out that 
working with the simpler set of equations helps reducing 
the risk of errors in the numerical implementation. 

In three dimensions and for B lmp = Fm = 0, but he- 
lical or non- helical kinetic forcing via .Fk, this system of 
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equations is capable of reproducing essential features of tur- 
bulent dynamos like initial exponential growth and subse- 
quent saturation; see, e.g., Brandenburg (2001) or Haugen 
et al. (2004). 

If .Bimp 7^ or Fm / we are no longer dealing with 
a dynamo problem in the strictest sense. A discussion of 
dynamo processes can still be meaningful if B{ mp = and 
the magnetic forcing does not give rise to a mean electro- 
motive force £q. A possibility to accomplish this is Fk = 
together with a magnetic forcing resulting in a Beltrami 
field bo, but any choice providing an isotropic background 
turbulence bo, Uo should be suited likewise. 

Note that the mean-field induction equation is still au- 
tonomous allowing for the solution B — 0. It then depends 
on properties of the background turbulence like chirality 
whether, e.g., the a effect or other mean-field effects render 
this solution unstable by enabling growing solutions. 

If we admit B lmp ^ 0, the two forcing terms can in 
principle be adjusted such that the resulting background 
turbulence is again isotropic and £o = 0. Should a grow- 
ing mean field then be observed, it can legitimately be at- 
tributed to an instability as B = B lmp is a solution of 
the mean field induction equation and the imposed field 
cannot grow. Thus both scenarios have the potential to ex- 
hibit mean-field dynamos although the original induction 
equation is inhomogeneous and the dynamo must not be 
considered as an instability of the completely non-magnetic 
state. 

Models of the latter type may well have astrophysi- 
cal relevance, because at high magnetic Reynolds numbers 
small-scale dynamo action is expected to be ubiquitous. 
Large-scale fields are still considered tobe the consequence 
of an instability, at least if there is no So or any other sort 
of "battery effect" . Magnetic forcing can be regarded as a 
modeling tool for providing a magnetic background turbu- 
lence when, e.g., in a DNS the conditions for small-scale 
dynamo action are not afforded. 

More generally, magnetic forcing and an imposed field 
provide excellent means of studying the a effect, the in- 
verse cascade of magnetic helicity, and flow properties in 
the magnetically dominated regime (see, e.g., Pouquet et 
al. 1976; Brandenburg et al. 2002; Brandenburg & Kapyla 
2007). 

3.2. Purely magnetic background turbulence 

Before taking on the most general situation of both mag- 
netic and velocity fluctuations in the background turbu- 
lence it seems instructive to look first at the case comple- 
mentary to that discussed in Section [2] That is, we assume, 
perhaps somewhat artificially, that the background velocity 
fluctuations vanish, i.e. Uo = 0, so that u = Ug. According 
to ([3]) we now find 



S = £g — Ug x b — u x b. (12) 

The modified momentum equation for the velocity fluctua- 
tions in a homogeneous medium reads (cf. (jlOp ) 

( )u 

— = J xb + j x B + T' + vV 2 u, (13) 



with T' = j B x b+j x b s -j s x b + j x b B . As (j x b )' 
needs to vanish in order to guarantee uq = 0, this could 



also be written as J 7 ' = j X b — j X 6. Unlike in the quasi- 
kinematic method there is now no longer any way to base 
a test-field method upon considering one of the fluctuating 
fields, here b, to be given (e.g. taken from a main run) while 
interpreting the other, here u, as a linear and homogeneous 
functional of the mean field. (This would work here, how- 
ever, in the second order correlation approximation, where 
J-' is set to zero.) 

3.3. General mean-field treatment 

The mean-field equations for U and B = curl A + B- imp 
obtained by averaging Eqs. ([TU]) and (fTTj) are 

ffFJ 

— = vV 2 U + JxB + F, (14) 
ot 

ffA _ 

— = nV 2 A + U x B + £, (15) 
ot 

where we have assumed that the mean forcing terms vanish. 
From now on we extend our considerations also onto the 
relation between the mean force J- = j x b and the mean 
field. In full analogy to the mean electromotive force we 
find this relation, considered as a functional of u, U and 
B, again to be linear and homogeneous in the latter and 
write, to start with, 

!F s = (t>B - ^VB. (16) 

Like ol and rj, the tensors and %p may depend on B. 
In the kinematic limit <p an d are expected to be non- 
vanishing only if bo 7^ 0. An analysis in SOCA, however, 
would also require Uo to get a non- vanishing result; 
see Appendix iBl Note that bo ^ allows T B to be linear 
in B, which would otherwise be quadratic to leading order. 
Consequently, the backreaction of the mean field onto the 
flow is no longer independent of its sign. 

As J-g is the divergence of the mean Maxwell ten- 
sor, it has to vanish in the homogeneous case, i.e. for 
homogeneous turbulence and a uniform field. Hence, for 
Equation (|16[) to be valid in physical space, <fi has then 
to vanish. However, in Fourier space we may retain relation 
(Til))) with lim fc _ > Q 4>(k) — (but not so for if)). On the other 
hand, in the physical space a description of employing 
the second derivative of B is likely to be more appropriate, 

i.e. 

T E = <f>* V 2 B - tpVB. (17) 

According to the expression for 4>{k), which is derived in 
Appendix [B] for Roberts forcing, Equation (|17p specified to 

would indeed be sufficient as long as there is sufficient 
scale separation between mean and fluctuating fields. In 
the following, we continue referring to <fi as introduced by 
Equation P^)l . 

The equations for the fluctuations are obtained by sub- 
tracting (fill) from (flU]), and (jTSj) from (jTTJ) , what leads to 

— = Jxb + jxB + F' + f K + vV 2 u, (18) 
= U xb + uxB + t + f M + r]\7 2 a, (19) 
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respectively, where !F' = j X b— j x b and £' = u x b it x b 
are terms that are quadratic in the correlations, while / k,m 
are just the fluctuating parts of the forcing functions. 

To arrive at a set of equations that are formally linear 
and allow for solutions as responses to a given mean field 
that are formally linear and homogeneous in the latter we 
make use of the split of all quantities into parts existing in 
the absence of B and parts vanishing with B introduced in 
Sect.[3J We write u = u +u B , a = a +a B and j = jo+Jg, 



For the parts vanishing with B we have 



T' s and & = £' 



£'^ and assume 

B 



that the forcing is independent of B. Equation ([T5]) and 
(fT^|) split consequently as follows (see Appendix [C] for an 
illustration) 



= vV 2 u +F + f K , 



dt 
dap 

at 



dan t — 

77V 2 a + [/x bo + e' + fu, 



du 



dt 
da s 



B = vW 2 u B + Jxb + jxB + F' s 



(20) 
(21) 

(22) 



dt 



iS -r)V 2 a B + U xb B + ux B + £' E . (23) 



Because of T'q = (jo x b )' and £' = (u x bo)', Eqs. ([20]) 
and (|2~T|) are completely closed. Furthermore, we have 

F' E = (jo xb B + j B xb + j B x b s )', (24) 
£' s = (uo x b B + u B x bo + u B x b s )' . (25) 

We can rewrite these expressions such that they become 
formally linear in u B and b B each in two different flavors: 

F'b = V xb B + 3b x b o)' = Uo x b B + j B x b)', (26) 
£' B = (u x b B + u B x bo)' = (u a x b B + u B x b)' . (27) 

Now we have achieved our goal of deriving a system of for- 
mally linear equations defining the part of the fluctuations 
that can be related to the mean field as response to the 
interaction with the given fluctuating fields it, Uq, b, and 
bo- 

Splitting mean force and electromotive force we find 



Fo = jo x b and £ = u x b 



(28) 



for the parts existing already with B = 0, due to a small- 
scale dynamo or magnetic forcing. Although it is hard to 
imagine that isotropic forcing alone is capable of enabling 
a non- vanishing £o, an additional vector influencing the 
otherwise isotropic turbulence may well act in this way. For 
example, using the second-order correlation approximation 
(SOCA) it was found that in the presence of a non-uniform 
mean flow U, with mean vorticity W = curl U, we have in 
ideal MHD (r) = u = 0) 



£o 



TU ■ 



Y u ° ■ 3° U + ~3~ u ° ■ b ° W 



(29) 



Here the index "0" refers to the fluctuating background 
uninfluenced by both the magnetic field and the mean 
flow. Beyond this specific result, too, one may expect that 
some cross correlation of the primary turbulences is crucial. 
(Yoshizawa 1990; Radler & Brandenburg 2010). 



Fb = 3 x b B + j B x b = jo x b B + j B x b, (30) 
£ B = u x b B + u B x bo = u x b B + u B x b. (31) 

We recall that for bo = (see Section [2]) , only the term 

u x bs = £^ occurs in the mean electromotive force and 
B b 

for u = (see Equation ([12"])) only u B x b = £^. For 

interpretation purposes, it is therefore convenient to define 

a correspondingly symmetrized version, 



Fb=3 *b s + j s xb- j B x b E 



£ B = u x b B + u B x b — u B x b B = £ B + £ B 



M | 7R 



with T B — —j B x b B and £j± = — u B x b B being resid- 
ual terms. Of course this split is only meaningful with a 
non- vanishing mean field in the main run. The correspond- 
ing transport coefficients might be split analogously. Note, 
however, that for an imposed field in, say, the x direction 
this is restricted to the (lj) components of the tensors. 



3.4. Test- field method 

In a next step we define the actual test equations starting 
from Eqs. @, ((Ml), © and ([27]). As tne y are already 
arranged to be formally linear when deliberately ignoring 
the relations between u B and u as well as between b B and 
b, respectively, we have nothing more to do than to identify 
B with a test field B T and b B ,u B with the corresponding 
test solutions b T ,tt T . Due to the ambiguity in (|26p and (12"T1) 
four different versions are obtained reading 



du 

da T 

dt 



= J T x b + j x B T + F' T + v\J 2 u T 



= U x b + ti x B 



7?V 2 a T , 



with 



(jxb T +j T xb )' 
or 

(joxb T +j T xb)\ 

(ux b T + u T x bo)' 
or 

(u Q x b T +u T x h)'. 



(32) 
(33) 

(34) 
(35) 



For mean force and electromotive force expressed by the 
test solutions we write correspondingly 



j x b T 4 


- j T x b 


oi- 




jo x b T 


+ j T x b, 


u x b T - 


- u T x bo 


or 





(36) 



(37) 



t Mo x b T + u T x b, 



and stipulate that the choice within Eqs. (|36l) and (13"T|) is 
always to correspond to the choice in Eqs. (|34[) and (|3"S"|) . 
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Table 1. Illustration of generating the four versions of 
the generalized test-field method by combining the different 
ways of representing T' T and £' T in Eqs. (|3~4"|) and (|3"5j) 





j x b T + j T x b 


jo x b T + j T x b 


u x b T + u T x bo 


ju 


bu 


u x b T + u T x b 


jb 


bb 



As we will make use of all four possible versions we label 
them in a unique way by the names of the fluctuating fields 
of the main run that enter the expressions for !F' T and £' T . 
Accordingly, we find by inspection of Eqs. (IMl) and (l3"5j) for 
the labels the combinations ju, jb, bu and bb; see Table Q] 

Now we conclude that for given u, 6, ito, and U the 
test solutions u T , b T are linear and homogeneousm the test 
fields B T and that the same holds for 'W T and £ T . Hence, 
the tensors a, tj, 4> and i/> derived from these quantities will 
not depend on the test fields, but exclusively reflect prop- 
erties of the given fluctuating fields and the mean velocity. 
If these are affected by a mean field in the main run the 
tensor components will show a dependence on B. Thus, like 
in the quasi-kinematic method, quenching behavior can be 
identified. We observe further that when using the mean 
field from the main run as one of the test fields, the corre- 
sponding test solutions b T and u T will coincide with b—bg 
and it — tto, respectively. 

Summing up, we may claim that the presented general- 
ized test-field method in either shape satisfies certain nec- 
essary conditions for the correctness of its results. But can 
we be confident, that these are sufficient, too? An obvi- 
ous complication lies in the fact, that in contrast to the 
quasi-kinematic method yielding the transport coefficients 
uniquely, we have now to deal with four different versions 
which need not be completely equivalent. Indeed we will 
demonstrate that the reformulation of the original prob- 
lem into Eqs. ([32]) and ([33]) with Eqs. ([34]) and ([31]) in- 
troduces occasional spurious instabilities. As we presently 
see no strict mathematical argument for the identity of the 
outcomes of all four versions, we resort to an empirical justi- 
fication of our approach. This is what the rest of this paper 
is devoted to. 



Remarks: (i) Applying the second order correlation approx- 
imation (SOCA) to the system (|32l) . ([33]) . that is, neglecting 
1F' T and £' T , melts the four versions down to one and thus 
removes all complications. 

(ii) In the limit B — > we have simultaneously u — > ito and 
b — > 6o, so again only one version remains. The method has 
then of course no longer any value for quenching considera- 
tions, but it still might be useful to overcome the limitations 
of SOCA. 

(iii) For b = the a T equation (|3"31 with the first version 
of (ES), i.e. 



£' T = (u x br)' 



(38) 



and correspondingly £t = u x b T , but (|32l) ignored, re- 
verts to the quasi-kinematic method. For comparison we 
will occasionally apply this method even when bo ^ and 
label the quantities calculated in this way with an upper 
index "QK". 



From now on we define mean fields by averaging over 
two directions, here over the x and y directions, that is, all 
mean quantities depend merely on z (if at all) and we get 
a ID mean-field dynamo problem. As a consequence, B z is 
constant and there are only two non- vanishing components 
of V.B, namely J x and J y so only the evolution of B x and 
By has to be considered. Moreover, £ z is without influence 
on the evolution of B. Hence, instead of Eqs. (|16l) and ([5]) 
we can write 



bijBj 



aijBj 



Vij J j i 



(39) 



where the original rank-three tensors xp and rj are degener- 
ated to rank-two ones. 

Only the four components of either tensor with i,j — 
1, 2 are of interest, thus altogether 16 components need to 
be determined. As one test field B T comprises two relevant 
components and yields an and an £ T , each again with 
two relevant components, we need to consider solutions of 
Eqs. (gU) and with Eqs. ([M]) and ([35]) for a set of four 
different test fields. 



Selection of test fields: The simplest choice are homoge- 
neous fields in the x and y directions, but they are only 
suited to determine the a. tensor. 

All four tensors can be extracted by use of test fields 
with either the x or the y component proportional to ei- 
ther cos kz or to sinfcz and the other vanishing (see, e.g., 
Brandenburg 2005b; Brandenburg et al. 2008a, 2008b; Sur 
et al. 2008). That is, B T is either Bf c = 6 tp cos kz or 
B PS = Sip sin/cz, p = 1, 2, where the superscript pq, q = c, s 
labels the test field and the subscript % refers to its com- 
ponents. By varying the wavenumber k, the wanted tensor 
components can in principle be determined as functions 
of k, but are then no longer allowed to be interpreted in 
the usual way (cf. Brandenburg et al. 2008a). Here we re- 
frain from doing so and fix k to the smallest possible value 
k = 2n/L z where L z is the extent of the computational 
domain in z direction. But even then we introduce some 
errors because the harmonic test fields do not belong to 
the class of mean fields for which the ansatzes Eqs. (ITBl 
and ([5]) are exhaustive (see Sect. [2]). We must be aware 
that the tensors calculated in this way are "polluted" by 
some contributions from terms with higher derivatives of 
B in £ and T '. To monitor these departures we compare 
the a and <fi tensors found with harmonic test fields with 
those obtained with uniform ones. 

For each pair of test fields (B pc , B ps ) we determine 2x4 
unknowns by solving the linear systems 



F\ = fajBj 



i ~rP q 



O i — Oil j -D j 



(40) 



q = c,s. Note that there is no coupling between the systems 
for p — 1 and p — 2. Inversion of the rotation matrix 



cos kz — sin kz 
sin kz cos kz 



(41) 



(with the angle kz) provides the solutions explicitly, hence 
we have 



t3 1 = R 



Here the superscript "t" indicates transposition. 



(42) 
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3.5. Forcing functions, computational domains and boundary 
conditions 

For testing purposes, a common and convenient choice is 
the Roberts flow forcing function, 

/ = akt + V x with * = cosk x x cos k y y, (43) 

and the effective forcing wavenumber ki = (k% + k 2 ) 1 / 2 . 
With the chosen averaging the Roberts forcing is isotropic 
in the xy plane. Furthermore, a is a parameter controlling 
the helicity of the flow: with a = it is non- helical while for 
a = 1 it reaches maximum helicity. If not declared other- 
wise, we will employ just maximally helical Roberts forcing. 

The Roberts forcing function can be employed for ki- 
netic as well as magnetic forcing, so we write /k,m = 
ATk.m/j where the Nk,m are amplitudes having the units 
of acceleration and velocity squared, respectively. Note that 
for a = 1 Eq. (|43p yields a Beltrami field, i.e., it has the 
property curl/ = k[f. Therefore, provided B lmp = 0, 
kinetic and magnetic forcing act completely uninfluenced 
from each other and create just a flow and a magnetic field 
having exactly the Roberts geometry. This is not the case 
for a ^ 1. 

The computational domain is a cuboid with quadratic 
base L x = L y = 2n while its z extent remains adjustable 
and depends on the smallest wavenumber in the z direction, 
ki z , to be considered. However, as the Roberts forcing func- 
tion is not z dependent, the runs from which only a and cj> 
are extracted were carried out in 2D with k\ z = 0. 

We choose here k x = k y = k% where ki = 1 is the 
smallest wavenumber that fits into x and y extent of the 
computational domain. For random forcing the domain is 
always cubic, i.e. L x = L y = L z = 2n. In all cases we 
assume periodic boundary conditions in all directions. The 
results presented below are based on revision r 13439 of the 
Pencil CodeQ, which is a modular high-order code (sixth 
order in space and third-order in time) for solving a large 
range of different partial differential equations. 

3.6. Control parameters and non-dimensionalization 

In cases with an imposed magnetic field, we set Bj mp = 
(B o ,0, 0). Along with it the forcing amplitudes Nk,m are 
the most relevant control parameters. The only remaining 
one is the magnetic Prandtl number, PrM = v jt]. 

It is convenient to measure length in units of the inverse 
minimal wavenumber kx, time in units of 1/rjk 2 , velocity 
in units of r/ki , and the magnetic field also in units of rjki . 
The forcing amplitudes -/Vk,m are given in units of r] kf and 
r) 2 kf, respectively. Results will also be presented in dimen- 
sionless form: onj and ipij m units of ijki, rjij in units of r), 
<pij in units of r\k\ if not declared otherwise. Dimensionless 
quantities are denoted by a tilde throughout. 

The intensities of the physically relevant actual and pri- 
mary turbulences are readily measured by the magnetic 
Reynolds number and the Lundquist number, 

Re M = u ims /r)k { , Lu = b rms /rjkf , (44) 

where it rms and 6 rms are the rms values of fluctuating veloc- 
ity and magnetic field, respectively, and fcf is the effective 
forcing wavenumber. 

1 http : //pencil-code . googlecode . com 



4. Results 

Throughout this section we set the mean flow U to zero. An 
important criterion for the correctness of the generalized 
test-field methods is the agreement of their results with 
those of the imposed-field method. In most cases we checked 
for this criterion. Of course, the imposed-field method is 
only applicable if the actual mean field in the main run is 
uniform. If this is not the case, we are in some instances still 
able to perform validation by comparing with analytical 
results. 

4.1. Zero mean magnetic field 

In this section we assume that the mean field is absent or 
weak enough as not to affect the fluctuating fields markedly, 
that is, u ss Mo, b sa b - In particular it can then not render 
the transport coefficients anisotropic. Therefore, we denote 
by a and rjt simply the average of the first two diagonal 
components of a and T], i.e. a = (an -I- a22)/2 and i] t = 
{Vn +%2)/2, respectively. If not specified otherwise we set 
-Bimp = 10~ 3 or zero. Furthermore, we take PrM = 1, i.e. 

V = T). 

4.1.1. Purely hydrodynamic forcing 

In order to make contact with known results, we consider 
first the case of the hydrodynamically driven Roberts flow. 
In two dimensions, no small-scale dynamo is possible, hence 
bo = and itorms = N^/vk 2 . In three dimensions, how- 
ever, this solution could be unstable, but we have not 
yet employed sufficiently large ReM allowing for that. For 
ReM -C 1, a is given by (Brandenburg et al. 2008a) 

a/a 0K = Re M /[l + (k z /k { ) 2 }, a aK = -u lms /2, (45) 

where k z is the wavenumber of the harmonic test fields. The 
minus sign in ao accounts for the fact that the Roberts flow 
has positive helicity, which results in a negative a. 

Making use of the quasi-kinematic method, as well as of 
all four versions of the generalized method, we calculated 
a for Nm = 0, k z = (2D case) and values of iVx rang- 
ing from 0.01 to 100 with a ratio of 10; 5 rms grows then 
from 0.005 to 50. Figure [T] shows a/ao versus ReM (solid 
line) . Here the data points for all methods are indistinguish- 
able. All results also agree with those of the imposed-field 
method. 

Agreement with the SOCA result ([4"5"j) (dotted line) ex- 
ists for ReM <§C 1. For ReM > 1, this is not applicable, be- 
cause dropping the £' T term is then no longer justified. The 
SOCA values are nevertheless numerically reproducible by 
the test-field methods when ignoring the T' T and £' T terms 
in Eqs. (l3"2"j) and (|3"3"|) : see the diamond-shaped data points 
in Figure [TJ 

Corrections to the result (1431) with the £' T term retained 
were computed analytically by Radler et al. (2002a, b). The 
corresponding values are again well reproduced by all fla- 
vors of the generalized test-field method as well as by the 
imposed-field method. 

In the first line of Table [U we repeat the result for 
A^k = 1 and added that for test fields with the wavenum- 
ber k z = 1, from where we also come to know the turbulent 
diffusivity r/ t . Note the difference between the a values for 
k z — 1 and k z — roughly given by the factor y2 from 
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Fig. 1. a/a K vs. Rcm for purely kinetic Roberts forc- 
ing with PrM = 1 and k z — (2D case) from the quasi- 
kinematic and all versions of the generalized method (solid 
line with squares). Note the full agreement with Eq. P5|) 
(dotted line) for ReM "C 1. Diamonds indicate the results 
of the test-field methods with the T' T and £' T terms in Eqs. 
(j3"2")) and ([55)) neglected, again coinciding with (|4"5)) . 



Fig. 2. a/aoM versus Lu for purely magnetic Roberts forc- 
ing with PrM = 1 and k z = (2D case) from all versions of 
the generalized method (solid line with squares). Note the 
full agreement with Equation (|4"6")) (dotted line) for Lu <C 1. 
Diamonds give the results of the test-field methods with the 
£' T and !F' T terms in Q38[) neglected, again coinciding with 



for k z = kt = 1. Additionally, the results of the quasi- 
kinematic method for k z = 1, a® K and r)^ K , are shown. As 
expected, they coincide completely with a and rjt- 

4.1.2. Purely magnetic forcing 

Next we consider the case of purely magnetic Roberts forc- 
ing i.e. ./Vk = 0. Due to its Beltrami property, bo cx /, 
Mo = is a solution of Eqs. (|2T)|) and ([2~Tj) . A bifurcation 
leading to solutions with Uq ^ cannot be ruled out gen- 
erally, but was never observed. Thus we have for the rms 
value of the magnetic vector potential ao rm s = Nyi/r]k^, 
hence &orms = Nyi/rjki. The appropriate parameter for ex- 
pressing the strength of the fluctuating field(s) is now the 
Lundquist number and the corresponding analytic result 
for Lu -C 1 reads 

a/a 0M = (Lu/Pr M )/[l + (k z /k t )\ a QM = 36 rms /4 (46) 

(for the derivation see Appendix [D]) . It turns out that the 
sign of a coincides now with that of the helicity of the forc- 
ing function. Again, we consider first the two-dimensional 
case with k z = 0; see Figure [5J In analogy to purely hy- 
drodynamic forcing we find for Lu <C 1 agreement between 
all versions of the generalized test-field method (solid line 
with squares) with Equation (|46l) (dotted line). For higher 
values the SOCA versions (see Sect. 14.1. lj) accomplish the 
same; see diamond data points. Note, that for the last data 
point with Lu = 7 it was necessary to lower the strength of 
the imposed field to Bi mp /i]ki — 10~ 4 , because otherwise 
the solution of the main run becomes unstable and changes 
to a new pattern. 

The second line of Table[5]repeats the result for Nm = 1 , 
again amended by those for k z = 1 and the results of the 
quasi-kinematic method which is obviously unable to pro- 
duce correct answers. This is because the mean electromo- 
tive force is now given by itg x bo. which is only taken into 
account in the generalized method. Note further that rjt is 
positive both for hydrodynamic and magnetic forcings. 



Table 2. Kinematic results for a and fy t for purely hydrody- 
namic (Am = 0) purely magnetic (Ak = 0), and hydromag- 
netic Roberts forcing, PrM = 1. The wavenumber of the test 
field is k z = 1, except in the third column where k z = 0. Its 
results agree with those of the imposed-field method. cfi K 
and ?7j refer to the quasi-kinematic method. 

jV K Nm &(k z = 0) a a QK ft ?? t QK 



4.1.3. Hydromagnetic forcing 

As already pointed out in Sect. 13.51 in the absence of a mean 
field, for simultaneous kinetic and magnetic Roberts forcing 
there is a solution of Eqs. (|2TJ)) and (|2"Tj) consisting just of 
the solutions Uq and 6o of the system forced purely hydro- 
dynamically and magnetically, respectively. Again, a bifur- 
cation leading to another type of solution cannot be ruled 
out, but was not observed. Only within SOCA, however, 
the decoupling of «o and 6o lets the value of a for hydro- 
magnetic forcing be purely additive in the values for purely 
hydrodynamic and purely magnetic forcings. We denote the 
two latter by = a(bo = 0) and a m = oc(uq = 0), re- 
spectively. When abandoning SOCA, the terms (lig x bo)' 
and (jo x bg + jgX bo)' in Equations ((22]) and (|23|) provide 
couplings between Ug and bg and give rise to an additional 
"magnetokinetic" part, of a defined as c* m k = a cx^ — a m . 
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-0.0081 


0.0664 


-0.4734 


0.6604 


0.3086 
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0.4715 
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0.9428 


0.4714 
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Note that we use here lower case subscripts k, m, and 
mk to distinguish from the split introduced at the end of 
Section r3.3l which applies to the nonlinear case. In contrast, 
the occurrence of a m k is a purely kinematic effect. While 
ctk and a m are, to leading order (and hence in SOCA), 
quadratic in the respective background fluctuations the 
magnctokinctic term is of leading fourth order and is rep- 
resentable in schematic form as a m k oc (uqxPq boGto)- 

Lines 5 and 8 of Table [5] show cases with hydromag- 
netic forcing and amplitudes adjusted such that we would 
have fik = a m = 1 if SOCA were valid. In either case the 
preceding two lines present the corresponding purely forced 
cases. Lines 6 to 8 refer to the SOCA version of the test- field 
methods. It can be clearly seen that the results are additive 
only in the latter case. The value of a m k as inferred from 
lines 3 to 5 is —0.1 resulting in a considerable reduction of a 
in comparison with the purely additive value. This is owed 
to the strong forcing amplitudes leaving the applicability 
ranges of SOCA far behind. 

Figure [3] shows a m k for equally strong velocity and mag- 
netic fluctuations as a function of ReM = Lu together with 
«k, «m, Ok + «m and the resulting total value a. Note 
the significant difference between the naive extrapolation 
of SOCA, a m + ak, and the true a. In its inset the figure 
shows the numerical values of a m k in comparison to the 
result of a fourth order calculation a m k = — V / 2/64u;? ms &r ms 
(for the derivation see Appendix [E|) . Clearly, the validity 
range of this expression extends beyond ReM = Lu = 1 
and hence further than the one of SOCA. It remains to be 
studied whether the magnetokinetic contribution has a sig- 
nificant effect also in the more general case when Uo |f &o • If 
so, considering a to be the sum of a kinetic and a magnetic 
part, as often done in quenching considerations, may turn 
out to be too simplistic. Likewise one may wonder whether 




0.1 1.0 



Fig. 3. a versus ReM = Lu for hydromagnetic Roberts 
forcing with PrM = 1 and k z = (2D case). Along with 
the total value the constituents ,ak, a m and a m k as well as 
CKk + oe m are shown. Note the sign change in a at ReM ~ 
0.54. The inset shows a m k in comparison to the result of a 
fourth order analytical calculation (solid line). 

closure approaches to the determination of transport coef- 
ficients supposed to be superior to SOCA can be successful 
at all as long as they do not take fourth order correlations 
into account properly. 



For the tensors <j> ar| d V'. which turn out to show up with 
simultaneous hydromagnetic and magnetic forcing only (in 
addition, <p requires z-dependent mean fields) we have of 
course again isotropy, (pu = 4>22 = <(>, ipu = ^22 = ip- 

As a peculiarity of the Roberts flow, ip vanishes in the 
range of validity of SOCA if the helicity is maximum (a = 1 
in (02]) )• For this case the first three panels of Figure S] show 
the dependences a(k z ), T]{k z ) and <j>{k z ) with different val- 
ues of uorms = &0rms (data points, dotted lines). The last 
panel shows ip{k z ) for o = 0.5 and the same forcing am- 
plitudes as before. As explained above, Uq and bo can now 
no longer be forced independently from each other. Hence 
both fields can not show exactly the geometry defined by 
(|43[) and u rms and 6 rms diverge increasingly with increasing 
forcing. 

As demonstrated in Appendix|Bj 4>(k z ) cx k z /(k z + kf), 
a{k z ),r](k z ),^(k z ) cx l/(k 2 z + kf) in the SOCA limit. For 
comparison these functions are depicted by solid lines. Note 
the clear deviations from SOCA for ReM = Lu = 5, partic- 
ularly in a. Note also that the expression for ip was derived 
under the assumption that the background has the geome- 
try (|4"3]) . It is therefore not applicable in a strict sense. The 
clear disagreement with the values of ip from the test-field 
method for the high values of ReM and Lu are hence not 
only due to violating the validity constraint ReM <C 1. 

4.2. Dependence on the mean field 

We now admit dynamically effective mean fields and hence 
have to deal with anisotropic fluctuating fields u and b 
which result in anisotropic tensors a, r], cp and ip. For the 
chosen forcing, B is the only source of anisotropy in the xy 
plane, so a. has to have the form 

ctij = aiSij + a 2 B i B : j, i,j = 1,2, 

with B the unit vector in the direction of B (here the x 
direction). We obtain then an = a± + a 2 and a 22 = Q-x- 
Of course, the tensors 77, cf> and tp are built analogously. 
Clearly, irrespective whether the forcing is pure or mixed 
the effects of Bi mp prevent the fluctuating u and b from 
having Roberts geometry . 

For vanishing magnetic background, 60 = 0, the gener- 
alized methods are still expected to give results coinciding 
with those of the quasi-kinematic one, but with bo 7^ we 
will leave safe mathematical grounds and enter empirical 
work. 

4.2.1. Purely hydrodynamic forcing 

In this case we have £g = u x bg = £^ and all flavors of 
the generalized method have to yield results which coincide 
with those of the quasi-kinematic method. This is valid to 
very high accuracy for the ju and bu versions and somewhat 
less perfectly so for the bb and jb versions. We emphasize 
that the presence of £?i mp , although being solely responsible 
for the occurrence of magnetic fluctuations, does not result 
in a failure of the quasi-kinematic method as one might 
conclude from the model used by Courvoisier et al. (2010). 

Figure [5] presents the constituents of a. as functions of 
the imposed field in the 2D case with PrM = 1 We may 
conclude from the data that a 2 is negative and approxi- 
mately equal to aff. For values of B lmv /i]ki > 5 its mod- 
ulus approaches a 22 = &i and thus gives rise to the strong 
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Fig. 4. a(k z ), rj t (k z ), and (j>(k z ) for hydromagnetic Roberts forcing with PrM = 1 and a — 1 (left three panels), likewise 
t/j(k z ), but for a = 0.5 (rightmost panel). Solid lines correspond to SOCA results, cf. Appendix[Bj Curve labels refer to 
ReM = Lu or (ReM,Lu). 
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Fig. 5. an (solid line, filled circles) and a?22 (dashed line, 
open triangles) as functions of the imposed field strength 
-Bimp, compared with —a^ (dotted line, small dots) and 
0:22— aii = — oil (dotted line, open circles) for purely kinetic 
Roberts forcing with 2Vk = 1 and PrM = 1- «i! ~ a2 
throughout. Note that «ok < and that the a symbols 
in the legend refer to the normalized (hence sign-inverted) 
quantities. 



quenching of the effective a — an. Indeed, a(B imp ) can be 
described by a power law with an exponent —4 for large 
-Bimp- By comparing with computations in which the non- 
SOCA term was neglected, we have checked that this dis- 
crepancy is not a consequence of SOCA. This is at odds 
with analytic results predicting either a oc B~ 2 (Field et 
al. 1999; Rogachevskii & Kleeorin 2000) or oc B~ 3 (Moffatt 
1972; Riidiger 1974). Sur et al. (2007) suggested that this 
difference was due to the fact that the flows are either time- 
dependent or steady. However, inspecting their Figure 2, 
their numerical values for a do exhibit the B~ A power law. 
They also found that a a M , defined similarly to our an, 
increases quadratically with B for weak fields and declines 
quadratically for strong fields (Sur et al. 2007). This is in 
turn in agreement with our present results. 



4.2.2. Purely magnetic forcing 

Here, the mean electromotive force is simply £ — = 

itjj x b. This is true as long as significant velocities in the 
main run occur only due to the presence of a mean field, 
that is, as long as uq = (see above). While B is weak, 
£ is approximately Ug x 5o- However, one could speculate 
that, if the imposed field reaches appreciable levels, i.e., 
if u is sufficiently stron g, £ ca n, with good accuracy, be 
approximated by ~£^ = u x bg. Since the quasi- kinematic 
method takes just this term into account, it should then 
produce useful results. 

In Figure [6] we show the rms values of the resulting 
magnetic and velocity fields as functions of the imposed 
field strength for Am = 1, corresponding to Lu = 1/2 if 
-Bimp = 0. The data points can be fitted by expressions of 
the form 



l+SLpAB*' 6 °™s l + B? mp /i?r 



(47) 



where B* 1.8 A^m- Note, that indeed the velocity fluc- 
tuations become dominant over the magnetic ones for 
-B imp /f?fci > 2. 

The resulting finding, as shown in Figure [71 is com- 
pletely analogous to the one of Sect. 14.2.11 but now we see 



-«2 approaching a22 = c*i with increasing B 



imp ■ 



Hence, the supposition that the quasi-kinematic method 
could give reasonable results for strong mean fields has not 
proven true as is not approaching an, despite the dom- 
ination of u rms over 6 rms . Instead, the values from the quasi- 
kinematic method have the wrong sign and deviate in their 
moduli by several orders of magnitude. 

In Table [3] we compare, for different values of B lmp , 
the values of an and 022, obtained using the generalized 
test-field method, with those of a^± and those from the 
quasi- kinematic method, ajf and a^ 2 K , where the entire 
dynamics of bg has been ignored. Note, again, that the re- 
sults of all four version of the generalized test-field method 
agree with each other. 
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Fig. 6. Root-mean-square values 6 rms (filled circles) and 
u rms (open circles) as functions of the imposed field strength 
Si m p for purely magnetic Roberts forcing Nm — 1 with 
Pi'm = 1. The solid and dashed lines represent the fits given 
by Equation (|4"7|t . 
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Fig. 8. an (solid line, filled circles) as function of the 
imposed field strength B imp , compared with — (dotted 
line, small dots), an (dash-dotted line, open circles) and 
a^i (dotted line, open squares) for hydromagnetic Roberts 
forcing with Am = Nk. = 1 and PrM = 1. The inset shows 
ct22 (dashed line, open triangles) compared to an. 
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Fig. 7. an (solid line, filled circles) and «22 (dashed line, 
open triangles) as functions of the imposed field strength 
Bimp, compared with — afi (dotted line, small dots) and 
a22 — cen = —a2 (dotted line, open circles) for purely mag- 
netic Roberts forcing with Nm = 1 and Pi'm = 1. Note that 
a^i f=a a2 throughout. 

4.3. Hydromagnetic forcing 

Given that the a effect can be sensitive to the value of 
Pi'm, we study an and 022 as functions of PrM, keeping 
Lu/ReM = 1 and B- lrnv /vki = 1 fixed. The result is shown 
in Figure [9] 

In the interval 2 < PrM < 10 the a coefficients ex- 
hibit a high sensitivity with respect to PrM changing even 
their sign at PrM ~ 0.7 and 2, respectively. Note also the 
occurrence of a minimum and the subsequent growth for 
Pr M > 4. 

Analogously to Figures [5] and [7] we show in Figure [H] the 
constituents of a versus B lmp . Note that we have used here 
chokRcmo + aoMLuo > for normalizing a. which is the 
kinematic value of an = aii for k z = and small itorms, 
Vms; see Equations |g5J) and (gBJ), Sect. 14.1.31 



Fig. 9. Dependence of an and 022 on PrM for hydromag- 
netic Roberts forcing with Lu/Pi'm = 1 and B lmp /vki = 1. 



It can be observed that oq\ at first dominates over —afi, 
but at Bimp/riki « 10 their relation reverses. Remarkably, 
the ratio of the moduli reaches, for high values of B lmpi just 



Table 3. Dependence of the diagonal components of the a 
tensor on B lmp for Nm = 1, PrM = 1 using the generalized 
method (an and (S22) together with the kinetic contribu- 
tion an ano - t ne results from the quasi-kinematic method 
{a n K and c$ 2 K ). 
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Fig. 10. Convergence of an the ju and jb versions of the 
generalized method with the result of the imposed-field 
method and exponential divergence of the versions bu and 
bb for Pr M = 1 with Nk — Nm — 1, B = 1, k z = and 
a resolution of either 32 2 mesh points (upper panel) or 64 2 
mesh points (lower panel). Note the significant improve- 
ment of the converging methods' agreement by doubling 
the resolution: The deviation is changing from ps 2.5% to 
ps 0.05%, that is, by a factor close to 2 6 suggested by the 
sixth order of the difference scheme. 



the inverse of that for low values. The strong quenching of 



-Q? x 



an is now a consequence of «5\ approaching 
In complete agreement with the former two cases with pure 
forcings, —an is proportional to B^ . However, we see a 
deviating behavior of a22(-Bimp) as it is no longer following 
a power law. 



4.4. Convergence 

In most of the cases the four different versions of the gener- 
alized method, (see Table [1]) give quite similar results. For 
purely hydrodynamic and purely magnetic forcing there is 
agreement to all significant digits. The agreement becomes 
somewhat less certain when there is hydromagnetic forc- 
ing, i.e. Nk ^ 0, Nm ^ 0. In general, however, agreement 
is improved by increasing the numerical resolution. 

Yet another complication arises when Bq ^ 0, because 
then some of the versions display exponentially growing 
solutions; see Figure ITOl We see no other explanation than 
that by the deliberate rearrangements leading to (|55"|) 
a potential for spurious instabilities was introduced. 



Their real occurrence, however, depends obviously on 
intricate properties of the fluctuating fields from the main 
run u and b. As the test equations are linear, unlimited 
exponential growth results (at least for stationary fields in 
the main runs). We suppose that if one could remove the 
unstable eigenvalues of the homogeneous version of the sys- 
tem (|32|) (133)) arbitrarily from its spectrum the solution of 
the inhomogeneous system would just be the correct one. 

5. Discussion 

The main purpose of the developed method consists in 
dealing with situations in which hydrodynamic and mag- 
netic fluctuations coexist. The quasi- kinematic method can 
only afford those constituents of the mean-field coefficients 
that are related solely to the hydrodynamic background Uq, 
but the new method is capable of delivering, in addition, 
those related to the magnetic background b . Moreover, it 
is able to detect mean-field effects that depend on cross 
correlations of Uo and bo. We have demonstrated this with 
the two fluctuations being forced externally to have the 
same Roberts-like geometry. With respect to a we observe a 
"magneto-kinetic" part that is, to leading order, quadratic 
in the magnetic Reynolds and Lundquist numbers and is 
capable of reducing the total a significantly in comparison 
with the sum of the a values resulting from purely hydro- 
dynamic and purely magnetic backgrounds. The tensors cf> 
and ip which give rise to the occurrence of mean forces 
proportional to B (or V 2 I?) and J are, to leading order, 
bilinear in ReM and Lu. In nature, however, external elec- 
tromotive forces imprinting finite cross-correlations of Uq 
and bo are rarely found. Therefore the question regarding 
the astrophysical relevance of these effects has to be posed. 
Given the high values of ReM in practically all cosmic bod- 
ies, small-scale dynamos are supposed to be ubiquitous and 
do indeed provide hydromagnetic background turbulence. 
But is it realistic to expect non-vanishing cross-correlations 
under these circumstances? 

Let us consider a number of similar, yet not completely 
identical turbulence cells arranged in a more or less regular 
pattern. We assume further that there is some asymme- 
try between upwellings and downdrafts such that, say, the 
downdrafts are more efficient in amplifying magnetic fields 
than the upwellings (Nordlund et al. 1992; Brandenburg 
et al. 1996). As dynamo fields are solutions of the homoge- 
neous induction equation and the Lorentz force is quadratic 
in B, bilinear cross-correlations, UQiboj, obtained by aver- 
aging over single cells can be expected to change their sign 
randomly from cell to cell provided the cellular dynamos 
have evolved independently from each other. Consequently, 
the average over many cells would approach zero and the 
aforementioned effects would not occur. In contrast, cross- 
correlations that are even functions of the components of bo 
and their derivatives, were not rendered zero due to polar- 
ity changes in the dynamo fields (e.g. the magneto-kinetic 
a). 

However, the assumption of independently acting cellu- 
lar dynamos can be put in question when the whole pro- 
cess beginning with the onset of the turbulence-creating 
instability (e.g. convection) is taken into account. During 
its early stages, i.e. for small magnetic Reynolds numbers, 
the flow is at first unable to allow for any dynamo action, 
but with growing amplitude the large-scale dynamo can be 
excited first to create a field that is coherent over many tur- 
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bulence cells. With further growth of its amplitude the (hy- 
drodynamic) turbulence eventually enters a stage in which 
small-scale dynamo action becomes possible. The seed fields 
for these dynamos are now prevailingly determined by the 
mean field and due to its spatial coherence the polarity of 
the small-scale field is not settling independently from cell 
to cell, thus potentially allowing for non-vanishing cross- 
correlations. 

But even if one wants to abstain from employing the in- 
fluence of a pre-existing mean field it has to be considered 
that neighboring cells are never exactly equal. Thus, in the 
course of the growing amplitude of the hydrodynamic back- 
ground in some of them, the small-scale dynamo will start 
working first hence setting the seed field for its immediate 
neighbors. It is well conceivable that the field polarity ini- 
tiated by one of the early starting cells "cascades" to more 
and more distant neighbors until this process is limited by 
the cascades originating from other early starting cells. The 
result could resemble the domains with uniform field orien- 
tation in ferromagnetic materials. Consequently, we arrive 
at a situation similar to the one discussed before, yet with 
less extended regions of coinciding field polarity. Instead of 
employing the idea of a pre-existing large scale dynamo in 
the cosmic object at hand one may even suppose that, given 
the smallness of the turbulence cells compared to the scale 
of the embedding surroundings, there is always a large-scale 
field, e.g. the galactic one, that is coherent across a large 
number of turbulence cells. 

In summary, cross-correlations and the mean-field ef- 
fects connected to them are to be considered realistic op- 
tions. Direct numerical simulations employing the scenar- 
ios discussed above should be performed in order to clarify 
the significance of these effects. This is equally valid for 
the cross-correlations related effects leading to £q (see Eq. 

In a recent paper, Courvoisier et al. (2010) discuss 
the range of applicability of the quasi-kinematic test-field 
method. Their model consists of the equations of incom- 
pressible magneto-hydrodynamics with purely hydrody- 
namic forcing. However, by imposing an additional uniform 
magnetic field B, together with the forced fluctuating veloc- 
ity a fluctuating magnetic field arises. It must be stressed 
that, following the further lines of their arguments, these 
fluctuations have to be considered as part of the background 
(uo,bo), that is, representing just those fluctuations that 
occur in the absence of the mean field. This follows from 
the fact that, when defining transport coefficients such as 
a, the field B is not regarded as part of the mean field 
(£?), in contrast to our treatment; see their section 2.(b). 
For simplicity they consider only the kinematic case and 
restrict the analysis to mean fields oc e lkzZ with k z — > 0. In 
their main conclusion, drawn under these conditions, they 
state that the quasi- kinematic test-field method which con- 
siders only the magnetic response to a mean magnetic field 
must fail for B ^ 0, that is bo ^ 0. We fully agree in this 
respect, but should point out that the method was never 
claimed to be applicable in that case; see Brandenburg et 
al. (2008c, Sect. 3), where it is mentioned that "as in al- 
most all supercritical runs a small-scale dynamo is opera- 
tive, our results which are derived under the assumption of 
its influence being negligible may contain a systematic er- 
ror." However, Courvoisier et al. (2010) overinterpret their 
finding in postulating that already the determination of 



quenched coefficients such as a(B) for bo = by means 
of the quasi-kinematic method leads to wrong results. The 
paper of Tilgner & Brandenburg (2008), quoted by them in 
this context, is just proving the correctness of the method, 
as do Brandenburg et al. (2008c). 

Our tensor ip is related to their newly introduced mean- 
field coefficient T by ipij = tkjzXi'ik- Unfortunately, an at- 
tempt to reproduce their results for T (and likewise for a) 
is not currently possible owing to our modified hydrody- 
namics. We postpone this task to a future paper. 

6. Conclusions 

Having been applied to situations with a magnetohydro- 
dynamic background where both Uq and bo have Roberts 
geometry, the proposed method has proven its potential for 
determining turbulent transport coefficients. In particular, 
effects connected with cross-correlations between Uq and bo 
could be identified and are in full agreement with analyti- 
cal predictions as far as available. No basic restrictions with 
respect to the magnetic Reynolds number or the strength 
of the mean field in the main run, which causes the nonlin- 
earity of the problem, are observed so far. As a next step, 
of course, the simplifications in the hydrodynamics we used 
will be dropped, thus allowing to produce more relevant re- 
sults and facilitating comparisons with work already done. 

Due to the fact that we have no strict mathematical 
proof for its correctness, there can be no full certainty about 
the general reliability of the method. As a hopeful indica- 
tion, in many cases, all four flavors of the method produce 
practically identical results, but occasionally some of them 
show, for unknown reasons, unstable behavior in the test 
solutions. Clearly, further exploration of the method's de- 
gree of reliance by including three-dimensional and time- 
dependent backgrounds is necessary. Homogeneity should 
be abandoned and backgrounds which come closer to real 
turbulence such as forced turbulence or turbulent convec- 
tion in a layer are to be taken into account. 

Thus, the utilized approach of establishing a test-field 
procedure in a situation where the governing equations 
are inherently nonlinear, although by virtue of the Lorentz 
force only, has proven to be promising. This fact encourages 
us to develop test-field methods for determining turbulent 
transport coefficients connected with similar nonlinearities 
in the momentum equation. An interesting target is the 
turbulent kinematic viscosity tensor, and especially its off- 
diagonal components that can give rise to a mean-field vor- 
ticity dynamo (Elperin et al. 2007; Kapyla et al. 2009), as 
well as the so-called anisotropic kinematic a effect (Frisch et 
al. 1987; Sulem et al. 1989; Brandenburg & von Rekowski 
2001; Courvoisier et al. 2010) and the A effect (Riidiger 
1980, 1982). Yet another example is given by the turbulent 
transport coefficients describing effective magnetic pressure 
and tension forces due to the quadratic dependence of the 
total Reynolds stress tensor on the mean magnetic field 
(e.g., Rogachevskii & Kleeorin 2007; Brandenburg et al. 
2010). 

Appendix A: Incompressibility 

The equations used in this paper had the advantage of sim- 
plifying the derivation of the generalized test-field method, 
but the resulting flows are not realistic because the pressure 
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and advective terms are absent. Here we drop these restric- 
tions and derive the test equations in the incompressible 
case with constant density. The full momentum and induc- 
tion equations take then the form 



dU 
dt 
OA 

~dt 



= U xW 
= U xB- 



hJxB + F K 
F M + rjV 2 A. 



vV 2 U 



VP; 



(A.l) 
(A.2) 



where W = curl U is the vorticity and P is the sum of gas 
and dynamical pressure and absorbs the constant density. 
The corresponding mean-field equations are 



dU 
dt_ 
dA 
~dt 



= U xW - 
= [7xBl 



J x B + F + vV 2 U - VP, 



r?V 2 A, 



(A.3) 
(A.4) 



where J : = uxw+jxb and £ = u x b, and the forcings 
were assumed to vanish on averaging. The equations for the 
fluctuations are consequently 



du 

at 

da 

at 



= Uxw + uxW + Jxb + jxB 
+ F' + F K + vV 2 u - Vp, 
= U x b + u xB + £' + F M + 



vV 2 a, 



(A.5) 
(A.6) 



where T' — (uxw + jx b)' and £' = (u x b)'. As above we 
split the fields and likewise the Equations ()A.5|) and (|A.6|) 



into two parts, i.e. we write u 
and arrive at 



U 



u B and a 



a 



duo 
dt 

da 
dt 



= U x w Q + u Q x W + F' + F K + iA7 2 « - Vp , 

(A.7) 



U xb + £' + F M +vV 2 ao, 



(A.8) 



and the equations for the B dependent parts 



a«R — 
-af =c/x ^ 



dan — 

— ^ = U xb R 

dt B 



xW + Jxb + jxB 



B 



V\J Ug — Vpg 

u x B c ' 



?]V 2 a E 



where T' = T' n 



F' s and £' = £' 



w Q +jo x bo)', £' = (u x b )', and 



£' s with T' 



F' s = (jo x b B + j B x b + j B x b B 

+ M X W B 



u B X W + u B X w B ) 



£'- 



(u xb B + u B xb 



u B x b B ) 



(A.9) 
(A.10) 

= (u x 



(A.ll) 
(A.12) 



We can rewrite these equations such that they become for- 
mally linear in u B and b B . Following the pattern utilized 
in Sect. 13.31 we find already for T' B four different ways of 
doing that. Together with the two variants in the case of 
£' B we finally obtain eight flavors of the test-field method 

where again in either case T and £ are to be constructed 



analogously to T' B and £' B . One of these flavors is defined 
by 



F' B = (u x w B + u B x w + j x b B + j B x b 



£' B = {ux b E 



x b )'. 



(A.13) 
(A.14) 



It is the one which comes closest to the quasi-kinematic 
test-field method, because there £' B — (u x b B )' . Next, we 

substitute B by a test field, B T , and u B and b B by the test 
solutions, u T and b T , i.e. 



du T — T 

dt 

da T „ 



x W 



dt 
where 

T' T = (u x w T 
£' T = (u x b T - 



+ F' T + vV 2 u T 
= TJ x b T + u x B T + £' T 



J T xb + j x B T 

Vp 1 , (A.15) 
(A.16) 



vV 2 a T , 



u 



xw +jxb T + j T x b )', (A.17) 



x b )'. 



(A.18) 



For the mean electromotive force and force the ansatzes ([5]) 
and (jTHJ) can be employed without change. Note, however, 
that the tensors <fi and ip now contain contributions from 
the Reynolds stress created by u B , that is eventually, by 
B. From the point of view of the tensorial structure of the 
relationship between B and T or £ the ansatzes © and 
(|16[) provide full generality as long as only the mean field 
and its first spatial derivative are to be included. That's 
why there are no separate terms with the mean velocity 
or its gradient tensor. Instead the latter play the role of 
problem parameters and all transport coefficients can of 
course be determined as functions of them. A separate task 
would consist in determining the tensors which appear in 
an analogous form of ((29)) in place of the scalar coefficients 
when a general anisotropic background is given. Then a test 
method with respect to U had to be tailored. 

Appendix B: Derivation of cf)(k z ), ip(k z ) 

Start with the stationary induction equation in SOCA 

r)V 2 b s + curl (u x B) = 0. (B.l) 

Assume u = uo^^f and_b = b Qlms f with / = f(x,y), 
curl/ = kff, f 2 = 1, B = Be lkzZ , and B XtV — const, 
B z = 0. Hence V 2 / = ~k 2 f. Then we can make the ansatz 
b B = b(x,y)e ikzZ with V 2 b = — k 2 b and get 

For the calculation of the mean force 



jo xb B + j B x b 



we need further 



j B = curl b B = curl (c ife ^b) = e ifc " z (curl b + \k z z x b) 



(B.2) 



r)(kf 



((£? • V)it + ik z (B ■ u )z) + ik z z x b^ 



(B.3) 
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Consequently 



T E = kfb x b E + j B x b 
1 1 



Tj k 2 



f 



k 2 



ik z [k { (B ■ u )z + kfii 0z B + z x (B ■ V)u ) x b 
+ k 2 z u 0z (z x B) x b 
1 1 



1,2 



(i*,*f(( 



ik z b 0z (B ■ V)m - £6 • (-B • V)« 
kl(u 0z b 0z B - zu 0z b ■ B) 



and with J — ik z z x B, that is, ik z Bk = Cfci3 Ji, = 1, 2, 

kf (ei3k£ml3Uo m bok — £ijk£kl3Uo z b j) J i 



m v k 2 + 



k' 2 



+ e/j3 -bo 



5i 3 b 



9«o \ -f 



9a;,- 



fc 2 ( u 0z b 0z Bi - d t3 u 0z boi B, 



Jl 



The tensors are hence 
1 k 2 



ll r\ k 2 



f 

1 1 



k 2 



(uo z b 0z 5a - u 0z b i8 i3 ,), 



rj k 



f 



k 2 



k{(u 0z b 0z S z i - uo z b ol Si 3 ) 



fef(l - S i3 )(u Qi b Q i - &u(u m b<Z\ + M02&02)) 

I ^ 3 



+ Qj3 yb 

<fe = 0*3 = 0. 



duo* 



bo ■ ^ — <J i3 



For fc 2 — > the tensor is proportional to fc 2 thus the corre- 
sponding mean force expressed in physical space by a convo- 
lution 4>o B, with 4> being the Fouricr-backtransformcd 0, 
can be approximated by a term oc d 2 B/dz 2 . For k z ^> fcf, 
however, the mean force can be represented by a term oc B. 
With Roberts geometry we have for a = 1 



1 

Tr{k 2 



"Orms borms , '0 = 0. 



All other components vanish, too. 

If, however, for the Roberts geometry < a < 1 the 
field / has indeed yet the property V 2 / = —k 2 f, but is no 
longer of Beltrami type. Instead, we have 

T 2 



curl / = <7kf(f 



1 



The tensor does not vanish any longer, but is now given 
by 



0H 

022 
012 



1 



fc 2 (l 



V(k 2 z + k 2 ) k f (l + <7 2 ) 

1 k 2 x (l-<J 2 ) 
V(k 2 + k 2 ) fc f (l + <7 2 ) 

021 = 0. 



^Orms^Orms i 



b 

^Ormsborms j 



Appendix C: Illustration of extracting a linear 
evolution equation from a nonlinear one 

To illustrate the procedure of extracting a linear evolution 
equation from a nonlinear problem, let us consider a sim- 
ple quadratic ordinary differential equation, y' = y 2 , where 
a prime denotes here differentiation. We split y into two 
parts, 2/ = 2/n +2/L; so we have 



y 



2/n 



2j/N2/L 



vl 



(C.l) 



In the last two terms we can replace ?/n+2/l by y, so we have 
2?/n2/l +2/l = (2/n + £/)i/L, which is now formally linear in j/l- 
Here, y corresponds to the solution of the 'main run'. Thus, 
at the expense of having to solve an additional nonlinear 
auxiliary equation, y^ = y 2 ^, we have extracted a linear 
evolution equation for t/l- Altogether we have 

,,2 



y = y , 

I 2 

Vn — y~N> 

vl = (2/n 



■u)ul, 



(C.2) 



where the last equation is linear in y\^. Note, that the system 
(|C.2[) is exactly equivalent to (jC.ip . i.e. no approximation 
has been made. 



with ^F'- dropped 



Appendix D: Derivation of ([351 

Consider the stationary version of 
(i.e. SOCA) and a uniform B, i.e., J = 

vV 2 u s + j x B = 0. (D.I 

Assume b = curia, diva = 0, hence j = — V 2 a. We get 

Ug = a x Bjv (D.2 

and further 

1 



(ug x b)i = —eu m eikjakb m Bj = onjBj 



that is, 



(a • bSij — aibj)/v. 



Isotropy results in 



a = an/3 = 2 a • bjiv. 

For b with Roberts geometry, however, we have a = an 
a 2 2 ^ a 33 , hence 



a = (a ■ b + a 3 b 3 )/2v = k f (a 2 + a\)j2v = 3b 2 . ms /Ak f u 
and with Lu = b rms /r]kf 



a = - 6 rms Lu/Pr M . 



(D.3) 



Adopt now a B depending on z only with B oc e lk * z , 
but 6, a still independent of z. Roberts geometry implies 
V 2 a = — k 2 a and V 2 ug = —(k 2 +k z )ug. Inserting in (jD.ll) 

(with the term proportional to J omitted) yields 



(fc 2 



k 2 )u s = k 2 f a x Bjv 



and comparison with (|D.2[) reveals that (|D.3I) has only to 
be modified by the factor l/(l + (k z /kf) 2 ). 
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Appendix E: Derivation of a m k in fourth order 
approximation 

We employ the iterative procedure described, e.g., in Radler 
& Rheinhardt (2007) to obtain those contributions to £g 
which are quadratic in norms an d &0rms an d expand for that 
purpose bg and itg into the series 





-<£>- 






U B 


= 4" 







with 

f] V 2 b t g ) = -curl(w x B) 
z/V 2 4 1} = -(jo x B + Jx b ) 



B 

72 



r,V z b^' = -curl(u xb| 1} + 1} x b ) 
^V 2 u« ^(jox&l;- 1 ' +4- 1) xb ), i = 2,... 



and 

oo oo 

% = E(-oX6f + «| ) x6o)=E^ ) 

i=l i=l 

In the following we assume B to be uniform and uq, bo to 
have Roberts geometry (|43p . The SOCA solutions b^ 1 -* and 

read 

From here on we switch to dimensionless quantities and set 
r) = v = 1, k x = k y = 1, fcf = \/2, = 1. So we have 

, (1) "Orms r . /X . i 

b^ = — - — [sini sm y, cosx cos y, — v 2 smi cosyj 
itg' = ^°^ ms [0 ; 2 cos x cos y , — y/2 sin % cos y] 



= K - w o rms [cos 2y, 0, sin 2y] + % s 
[cos 2y (cos 2a; + 2), sin2ysin2x, \/2 sin 2 j/ (cos 2x + 3) 

(3) (3) 

For b- and we present here only those parts which 

B 3 

eventually contribute to a m ^: 

■ (3) u 0rms^0rms r • A /n • 1 

b- = [sm x siny, cos x cosy, — 4v 2 sm a; cos y\ 



(3) "OrmsVms rn V% . i 

= ^ [0, cos x cos y, — — sm x cos yj 

Finally, 

7F(3) 



£•'■" ,(.•(! (3) , 2 1.2 

= M X b' s ' +U K g J X 6 = -W 0rms 6 rm S -^- 



i.e. 



V2 



2 ,2 v - 

u 0rms°0rms - ■ 



G4 



"("3") 

Note, that the contributions omitted in S- provide fourth 
order corrections to «k and a m . They result in dependences 
on Rcm and Lu that are weaker than the parabolic SOCA 
ones; see Fig. [31 
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